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Comment on "Effects of Point Defects on the 
Phase Diagram of Vortex States in High-T c Su- 
perconductors in the B || c Axis" 

Recently Nonomura and Hu (N-H) presented [1] simu- 
lations of a bond disordered three dimensional uniformly 
frustrated XY model, as a model for a point disordered 
type-II superconductor, using very strongly anisotropic 
couplings Jc/Jab — (A a fc/A c ) 2 = 1/400. As part of their 
phase diagram, they report a vortex slush (VS) phase 
that lies at disorder strengths e above the elastically dis- 
torted vortex line lattice (the Bragg glass (BG)) and is 
separated from it by a first order phase transition e c i(T). 
Increasing either disorder strength or temperature, the 
VS has another first order transition e C 2(T) to the vor- 
tex liquid (VL). The VS has neither superconducting 
phase coherence (longitudinal helicity modulus vanishes, 
T c = 0) nor translation order (liquid-like structure func- 
tion 5*(k)), and is distinguished from the VL by a sharp 
drop in the density of dislocations in the ab plane. Here 
we clarify the nature of the VS found by N-H, and argue 
that it is the result of an unphysical finite size effect. 

We have carried out simulations of N-H's model, using 
all the same parameters. We first cool to T = 0.08 at 
high e = 0.14, staying within the VL. Then we slowly 
decrease e to enter first the VS and then finally the BG. 
Then we slowly increase e back to 0.14. In Fig. la we 
plot our results for the helicity modulus, T c . In Fig. lb 
we plot S(K), where K is a reciprocal lattice vector of 
the vortex line lattice we find in the BG. Decreasing e, 
we find that T c rises from zero, indicating the onset of 
superconducting phase coherence, and S'(K) saturates to 
its low e value, indicating formation of an ordered vortex 
line lattice, both at e c \ ~ 0.06. Increasing e, however, 
we find that both phase coherence and the vortex line 
lattice persist to the higher value e C 2 ~ 0.10. Our values 
of e c i and e C 2 are reasonably close (given likely sample 
to sample fluctuations) to the values reported by N-H 
for the boundaries of the VS at T = 0.08. We therefore 
identify e c i < e < e C 2 as the VS and conclude that it is a 
region with considerable hysteresis and metastability. 
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when decreasing e. The different symbols represent dif- 
ferent values of K, corresponding to vortex lattices of 
different orientation. We see that in the VS, vortices 
in most planes have ordered into a 2D vortex lattice, 
however the orientation of this lattice varies in different 
layers. Decreasing e, the thickness of aligned layers in- 
creases, until below e c i in the BG, a single orientation 
fills the system. 




FIG. 1. a) Helicity modulus T c and b) structure function 
peak height S(K) vs. disorder strength e at T = 0.08. Solid 
(open) symbols are for decreasing (increasing) e. 



FIG. 2. Vortex structure function S'(K) computed layer by 
layer at heights z, for several different e, upon decreasing e. 
The different symbols correspond to different values of K, 
representing different vortex lattice orientations. 

For a particular random realization and system size L, 
it is possible that the above VS is indeed more stable than 
an ordered vortex line lattice for a region of e. However 
this cannot remain true as L increases. For an ordered 
vortex lattice of size L and thickness L z ~ L, the most 
energy that can be gained by adjusting to a particular 
set of random point pins scales as eJ a bL 3 ^ 2 . The energy 
loss by having a mismatch of vortex lattices between two 
adjacent planes will scale as J Z L 2 . For small J z , fixed L, 
such an adjustment may be favored. But as L increases, 
the energy cost always wins out and the vortex lattices 
will align for all layers. The vortex slush region will dis- 
appear, being replaced by a vortex line lattice, with a 
single melting transition to the liquid. 
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To get a clearer picture of the VS, we plot in Fig. 2 
5(K), computed layer by layer for the ab planes at differ- 
ent heights z, for several different values of e as obtained 
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